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Abstract. Using the method which retains the rotation symmetry of spin components in the paramagnetic 
state and has no preset magnetic ordering, spectral and magnetic properties of the two-dimensional t- 
J model in the normal state are investigated for the ranges of hole concentrations < x < 0.16 and 
temperatures O.Olt < T < 0.2t. The used hopping t and exchange J parameters of the model correspond 
to hole-doped cuprates. The obtained solutions are homogeneous which indicates that stripes and other 
types of phase separation are not connected with the strong electron correlations described by the model. 
A series of nearly equidistant maxima in the hole spectral function calculated for low T and x is connected 
with hole vibrations in the region of the perturbed short-range antiferromagnetic order. The hole spectrum 
has a pseudogap in the vicinity of (0,7r) and (n,0). For x « 0.05 the shape of the hole Fermi surface is 
transformed from four small ellipses around (±7r/2, ±7r/2) to two large rhombuses centered at (0,0) and 
(tt, 7t). The calculated temperature and concentration dependencies of the spin correlation length and the 
magnetic susceptibility are close to those observed in cuprate perovskites. These results offer explanations 
for the observed scaling of the static uniform susceptibility and for the changes in the spin-lattice relaxation 
and spin-echo decay rates in terms of the temperature and doping variations in the spin excitation spectrum 
of the model. 

PACS. 71.10.Fd Lattice fermion models - 74. 25. Ha Magnetic properties - 74.25.Jb Electronic structure 

1 Introduction the method takes proper account of the kinematic inter- 

action. The method is based on Mori's projection opera- 
The two-dimensional t-J model was proposed by Ander- tor technique [E] which allows one to represent Green's 
son Q] for the description of strong electron correlations f unc ti ns in the form of continued fractions and gives a 



in CuO-2 planes of perovskite high-T c superconductors. In way for ca l cu l a ting their elements. The residual term of 
Ref. the similarity of the low-energy part of its spec- the fraction is approximated by the decoupling which re- 
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O ■ trum with the spectrum of the realistic three-band Hub- duces t his many-particle Green's function to a product 

; bard model was demonstrated. Nowadays the t-J model f t he sought-for functions. Using the idea of Ref. [E] 

is one of the most frequently used models for the interpre- the decoupling is corrected by introducing a vertex cor- 

tation of experimental results in cuprates (for a review, rec tion, a multiplier which is determined from the sum 

see Ref. 0). Different numerical and analytical methods ru i e f t he considered problem, the constraint of zero site 

were used for the investigation of the model. Among these magnetization in the paramagnetic state. Test calculations 

methods are the exact diagonalization of small clusters @J with this method for small clusters and high temperatures 

E|, Monte Carlo simulations ©, density matrix renormal- demonstrated good agreement of the obtained results with 

ization group calculations 0, spin- wave [110 and mean- the exact diagonalization and Monte Carlo data. Notice 

field slave-boson approximations EDI In spite of the con- t hat close approaches which use equations of motion and 

siderable progress made towards the understanding of the Tserkovnikov's formalism were developed in Ref. Q3| to 

properties of the model, the basic issues of its behavior at consider magnetic properties of the 2D t-J model, 
moderate doping have not yet been completely resolved. 

Aiming at a description for this range of concentra- This article contains results of calculations carried out 

tion in Ref. JT] a new analytical method was developed with the method of Ref. m for a 20 x 20 lattice and pa- 

which has merits of retaining the rotation symmetry of rameters J = 0.1 eV, t = 0.5 eV where J and t are the 

spin components in the paramagnetic state and of the ab- exchange and hopping constants of the t-J model. Test cal- 

sence of any preset magnetic ordering. Since the descrip- culations with lattices of other sizes showed that the con- 

tion is carried out in terms of the Hubbard operators, sidered lattice is large enough to avoid finite-size effects in 
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spectral functions. The parameters chosen correspond to 
hole-doped cuprate perovskites ^3 • The considered range 
of hole concentrations < x < 0.16 spans the cases from 
light to moderate doping. The temperature range of the 
calculations O.Oli « 58 K < T < 0.2t w 1200 K embraces 
the range of temperatures used in experiment. 

The calculated Green's functions reveal a number of 
peculiarities in the hole and spin excitation spectrum. For 
low x and T, besides the spin-polaron peak [3], a series of 
less intensive and broader maxima is observed at higher 
frequencies in the hole spectral function. These maxima 
are nearly equidistant which allows us to connect them 
with vibrations of a hole in a region of the perturbed short- 
range antifcrromagnetic order. This region arises due to 
the hole movement. The lowest of these maxima with a 
dispersion which mimics the dispersion of the spin-polaron 
peak is apparently observed in lightly doped Ca2CuC>2Cl2. 
For low hole concentrations only the spin-polaron band 
crosses the Fermi level which leads to the Fermi surface 
consisting of four ellipses around (±7r/2, ±7t/2). For x « 
0.05 the second band, which arises below the Fermi level 
in the used hole picture, crosses the Fermi level and the 
Fermi surface acquires new elements - large rhombuses 
around (0, 0) and (rr, ir). However, due to the peculiar dis- 
persion of the spin-polaron band and its larger spectral 
intensity parts of this large Fermi surface near (0, tt) and 
(71", 0) are hidden which looks like a pseudogap, a decrease 
of the low-frequency spectral intensity in these regions of 
the Brillouin zone. The magnitude of the pseudogap - the 
energy distance between the spin-polaron peak and the 
Fermi level - decreases with x and at x « 0.12 the pseu- 
dogap disappears. The symmetry, magnitude and concen- 
tration dependence of the calculated pseudogap are similar 
to those observed in Bi-based cuprates. 

As mentioned, the method used has no preset magnetic 
ordering which opens the way to investigate whether the 
strong electron correlations described by the t-J model are 
responsible for the charge and magnetic inhomogeneities 
observed in some cuprates [IB] . In the considered ranges 
of hole concentrations and temperatures only the homo- 
geneous solutions were found which indicates that other 
interactions have to be invoked to explain these inhomo- 
geneities. 

The calculated spectrum of spin excitations contains a 
gap near (n, tt) which is directly connected with the corre- 
lation length of the short-range antiferromagnetic order. 
The dependence of the correlation length on x reproduces 
the relation observed experimentally in I^-zSr^CuO^ 
|17| . With increasing x the branch of spin excitations is 
destroyed near (0, 0), whereas at the periphery of the Bril- 
louin zone the excitations remain well defined. The calcu- 
lated magnetic susceptibility and its dependencies on x 
and T are similar to those observed in neutron scattering 
and NMR experiments. This similarity allows us to offer 
explanations for the observed scaling of the static uniform 
susceptibility and for the changes in the spin-lattice relax- 
ation and spin-echo decay rates in terms of the tempera- 
ture and doping variations in the spin excitation spectrum 
of the t-J model. 



For convenience main formulas of Ref . jllj , which were 
used in the calculations, are reproduced in Sec. II. In 
Sec. Ill and IV peculiarities of the hole and spin exci- 
tation spectra are discussed. The spin susceptibility, spin 
correlations, spin-lattice relaxation and spin-echo decay 
rates are considered in Sec. V. Our conclusions are given 
in Sec. VI. 



2 Main formulas 

The Hamiltonian of the 2D t-J model reads |3j 

H = ^ ' t Iim a\ la a ma -\- — ^ ' Jnm ( s n S m S n S m ) ' (-0 
rimer nm 

where a na — |ncr)(n0| is the hole annihilation operator, 
n and m label sites of the square lattice, a = ±1 is the 
spin projection, \na) and |n0) are site states correspond- 
ing to the absence and presence of a hole on the site. 
These states may be considered as linear combinations of 
the products of the 3d x 2_ y 2 copper and 2p a oxygen or- 
bitals of the extended Hubbard model |%lll8| . In this work 
we take into account nearest neighbor interactions only, 
t nm = -t J2 a <5n,m+a and J nm = J J] a ^n,m+a where the 
four vectors a connect nearest neighbor sites. The spin- 
\ operators can be written as s n = \ J2<r a\na) (na\ and 
s n = \wr){n,-a\. 

Due to the complicated commutation relations the di- 
agram technique for the operators et nCT , s n , and s n is very 
intricate In this case the use of Mori's projection 
operator technique |12) for the derivation of self-energy 
equations for Green's functions constructed from such op- 
erators is especially fruitful. In this way the following 
equations were found for the hole G(ki) = ((a kcr |a kcr )) = 

-i0{t)({a ka {t) )a tj and spin D(kt) = -#(t)K(t),*i k ]) 
Green's functions |11| : 

nflr s _ [4Ja(A + 1 + 7k )]- 1 77(ku;) + 4Jd( 7k - 1) 

(2) 

G(kw) = <j>[u) - £ k + n - U{kuj)}-\ 
where D(kcj) = exp(iu>t)D(\st)dt, 

akcr = N~ 1/2 exp(~^kn)a ncr , 

n 

4 = JV-V2^expHkn)4, 

n 

ako-(i) = exp(i7it)a k(T exp(— iTCt), N is the number of sites, 
H = H — ^ n I„, n is the chemical potential, X n = 
|n0)(n0|, the angular brackets denote averaging over the 
grand canonical ensemble, 7k = | J2a exp(ika), <f> = i(l + 
x), and 

col = 16J 2 a|C 1 |(l- 7 k)(^ + l + 7k), 

(3) 

£k = -(4# + 6Ci4>-H + SF^- 1 J) 7k . 
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In the present calculations the parameter of vertex cor- 
rection a, which improves the decoupling in the resid- 
ual term, is set equal to its value in the undoped case, 
a = 1.802 [this value differs slightly from that reported 
in Refs. [111120) due to an artificial broadening introduced 
in D(kuj), see below]. The parameter A which describes 
a gap in the spectrum of spin excitations at (tt, 7r) [see 
Eq. 0] is determined by the constraint of zero site mag- 
netization (sf) = in the paramagnetic state. The con- 
straint can be written in the form 



N 



E 



dwcothf — — ] BCkuj), 
\2TJ ' 



(4) 



where B(\tLu) = — 7r _1 Im D(\uu) is the spin spectral func- 
tion. Notice that in accord with the Mermin- Wagner the- 
orem in the considered 2D system the long-range an- 
tiferromagnetic ordering is destroyed at any nonzero T 
and, as will be seen below, for T = and x > x c « 0.02. 
The value of x and the nearest neighbor correlations C\ = 
( s l~ ls l+a) ana - = ( a i a i+a) are connected with Green's 
functions J3J) by the following relations: 

1 f°° 
x = — / dujnp(uj)A(k.Lu), 

1 f°° 
Fi = — ^Tk / dun F (w)AQsu), (5) 

C i = T^E W ^coth(— ) B(kw), 

k Ja 

where A(ku) = — 7r _1 Im G(kw) is the hole spectral func- 
tions and uf(<jj) — [exp(w/T) + 
The self-energies in Eq. (|2J) read 



ImTT(kw) 



16?rf 2 J 



(Z\ + l+ 7k ) V( 7k - 7k+k 



N 



k' 



dui'[np{uj + a/) — tif(u/)] 



Im E(kw) 



x A(k + k',uj + uj')A(\s'uj'), 
Wirt 2 



N(j> 



k' 



sgn(w')(7i 



k-k' 



'7k 



7k-k' + 7k 



(6) 



f 



7k' 



1 - 7k' 

x [n B (-Lo') + n F (uj - oj')] 
x A(k-k',w-w / )-B(kV), 

where ns(ti>) = [cxp(tu/T) — I] . The source of damping 
of spin excitations described by Eq. © is the decay into 
two fermions. Another source of damping, multiple spin 
excitation scattering, is considered phenomcnologically by 
adding the small artificial broadening — 2ryw k , f] = 0.02i to 
ImTI(kcj). The broadening —77 is also added to lmS(kuj) 
to widen narrow lines and to stabilize the iteration proce- 
dure. 



The same derivation for the transversal spin Green's 
function gives ((s^ 1 ^^ 1 )) = 2D(kt) indicating that the 
used approach retains properly the rotation symmetry of 
spin components in the paramagnetic state. 

For low hole concentrations and temperatures the band- 
width of the dispersion e k , Eq. is approximately equal 
to t which is much smaller than 8i, the bandwidth of un- 
corrected electrons. The reason for this band narrowing 
is the antifcrromagnetic alignment of spins. In this case 
the hole movement is accompanied by the spin flipping. 
In these conditions the hole dispersion is determined by 
the self-energy Z'(kcj) which is responsible for a consider- 
able energy gain for states in the spin-polaron band. For 
the parameters used this gain is approximately equal to 
2t. 



3 The hole spectrum 

Equations @~© form a closed set which was solved by 
iteration. Examples of the hole spectral functions for the 
cases of low and moderate doping are given in Figs. ^ and 
121 The frequency is measured from the chemical potential. 
For low x the shapes of the spectra are close to those ob- 
tained in the spin- wave approximation [Hj. As seen from 
Fig. n a t low doping the hole spectral function contains 
a series of nearly equidistant maxima. The first member 
of this series is the narrow peak with the highest inten- 
sity which is frequently termed the spin-polaron peak (in 
some of the spectra in Figs.n an( i|21 0m y its foot is shown). 
The dispersion of the maxima in the above spectra along 
the symmetry directions is shown in Fig. [21 In the spin- 
polaron band the lowest energy is reached at (tt/2,tt/2). 
The nearly equidistant location of the maxima is clearly 
seen in Fig. |3K- Besides, we notice that the second maxi- 
mum in the series has a dispersion which is close in shape 
to the dispersion of the spin-polaron peak. As follows from 
Figs.ElandEb, this maximum is retained at moderate dop- 
ing in some region of the Brillouin zone. 

The nearly equidistant position of the maxima in the 
low-concentration spectra allows us to connect them with 
vibronic states of a hole in the region of the perturbed 
short-range antifcrromagnetic order. The moving hole lea- 
ves behind it a trace of overturned spins. Such disturbance 
requires energy which in its turn leads to a restoring force 
acting on the hole and giving rise to its vibrations. This 
notion is close to the string picture developed in Ref. |2*2*] 
for the Ising model and used for the interpretation of the 
fine spectral structure in the t-J model on small lattices 

The two lowest maxima of the above-mentioned se- 
ries are apparently observed in lightly doped Ca2Cu02Cl2 
|2.'i| . The narrow and intensive spin-polaron peak is re- 
sponsible for the photoemission maximum with the low- 
est binding energy, while the second member of the series 
corresponds to the spectral maximum which is observed 
at approximately 600 meV higher binding energy. This 
energy difference is close to that shown in Fig. HJi. Be- 
sides, as noted in Ref. [221 , the dispersion of the second 
maximum mimics the dispersion of the spin-polaron peak, 
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-2 6 " 2 4 6 <a" 

ro/f Fig. 2. The hole spectral function for x = 0.12 and T = 0.02t. 

Fig. 1. The hole spectral function A(ku>) along the symmetry 
lines for x = 0.012, T = 0.02t and J/t = 0.2. The respective 
values of the wave vectors are indicated near the curves. Here 
and below the hole picture is used. The vertical dotted line 
indicates the position of the chemical potential. 
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Fig. 3. The dispersion of the maxima in the hole spectral 
function for T = 0.02t, x = 0.012 (a) and x = 0.12 (b). Points 
Y, M and 5* correspond to k = (0,7r), (7T,7r) and (n /2,n/2), 
respectively. 



which conforms with our calculations. However, it should 
be noted that our calculated dispersion of the spin-polaron 
band differs somewhat from that observed experimentally 
in lightly doped cuprates 23 , 24 . Although the shapes and 
magnitudes of the experimental and calculated dispersions 
along the nodal direction [from (0, 0) to (ir, it)] are close, 
along the boundary of the magnetic Brillouin zone [from 
(0, 7r) to (it, 0)] the magnitude of the calculated dispersion 
is an order of magnitude smaller than the experimental 
one. A possible reason for this discrepancy is the oversim- 
plified hole hopping term in Hamiltonian Q which takes 
into account the transfer between nearest neighbor sites 
only 1231125] . 

The spectral intensity below the chemical potential in 
Figs. and |3K is noteworthy. Already at low x in some 
region of the Brillouin zone this spectral feature looks 
like a weak maximum. It is this spectral intensity which 
provides a finite hole concentration until the spin-polaron 
band crosses the Fermi level at x « 0.04. 

With increasing x and T the high-frequency maxima in 
the nearly equidistant series are smeared, as seen in Fig.UJ 
Inspecting this figure one can see that for moderate hole 
concentrations two maxima - the spin-polaron peak and 
a less intensive and broader maximum - cross the Fermi 



(0,71) (7C,Jt) 



x=0 .043 




(0,0) (71,0) 




Fig. 4. The Fermi surface in the first quadrant of the Brillouin 
zone for x = 0.043, 0.08, and 0.12. T = 0.02t. 

level. This latter maximum originates from the mentioned 
spectral feature which appears below the chemical poten- 
tial at low x. The dispersions of the two maxima and some 
other spectral peculiarities are shown in Fig. Eh- The max- 
ima can be resolved in the spectrum only in some regions 
of the Brillouin zone. This is the reason why the curves in 
Fig. [3Jd terminate at some points of the symmetry lines. 

The evolution of the Fermi surface with doping is shown 
in Fig. 2] No Fermi surface exists for x < 0.04, since 
no maximum which can be identified with a quasipar- 
ticle excitation crosses the Fermi level (see Figs. ^ and 
As indicated above, for such x the occupied hole 
states are located at some distance from the level. At 
x w 0.04 the spin-polaron band crosses the Fermi level 
near (±7r/2, ±7r/2) which produces a Fermi surface con- 
sisting of four strongly elongated ellipses around these 
points. In Fig.0]the transverse size of one of these ellipses 
is much smaller than the used momentum step and there- 
fore the ellipse is depicted as a line segment. At x w 0.05 
the second band crosses the Fermi level. With this cross- 
ing the Fermi surface acquires new elements - in addition 
to the ellipses of the spin-polaron band there appear two 
large rhombuses centered at (0, 0) and (ir, ir) (see Fig. 01 
the case x = 0.08). One of these rhombuses can be con- 
sidered as a shadow image of the other, which arises due 
to the antiferromagnetic short-range order. Notice, how- 
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ever, that the rhombuses are slightly different in size. In 
such a manner the small Fermi surface which exists for 
x < 0.05 is transformed to the large surface for larger hole 
concentrations. With a further increase of x the Fermi 
level moves up from the bottom of the spin-polaron band, 
the size of the ellipses grows and their shape is somewhat 
changed. At x w 0.12 they reach the boundaries of the 
Brillouin zone and the parts of the Fermi surface con- 
nected with the spin-polaron band are transformed into 
two rhombuses the shapes and sizes of which are close to 
those connected with the second band. Thus, for x > 0.12 
the Fermi surface consists of two rhombuses centered at 
(0,0) and (tt, tt) (see Fig.EJ). 

The compound Bi 2 Sr 2 CaCu208+<5 is most intensively 
investigated in angle-resolved photoemission experiments. 
However, even for this crystal the Fermi surface topology 
is a topic of an intense controversy (see, e.g., Ref. |26j and 
references therein). The main reasons for this are broad 
spectral features in underdoped normal-state crystals and 
the weak dispersion around (tt, 0) (the so called flat bands 
or extended van Hove singularities, see Fig.Qb). We notice 
also that for moderate doping the spin-polaron band in 
all the area between the rhombuses lies less than 20 meV 
below the Fermi level. This energy is on the verge of ac- 
curacy of photoemission experiments. Nevertheless in the 
variety of experimental results on the Fermi surface of 
normal-state cuprates it is possible to find some features 
which resemble those shown in Fig. 01 The shadow Fermi 
surfaces (together with some other replicas of the main 
Fermi surface) are observed in Bi 2 Sr2CaCu208+<5 [2HI23- 
In some experimental conditions in this crystal and in 
La2- a; Sr a ;Cu04 the shapes of experimental Fermi surfaces 
|2lHl2"5] are close to that shown in Fig. H for x = 0.12. 
Notice however that in Ref. [23 such shape of the Fermi 
surface in I^-rrSr^CuC^ is ascribed to dynamic stripes 
presumed in this crystal. 

As indicated, in the range 0.05 < x < 0.12 the Fermi 
surface contains parts arising due to the crossings of the 
Fermi level by two different bands. The spectral maxima 
corresponding to these bands differ essentially in their 
intensity. In comparison with the spin-polaron peak the 
spectral maximum of the second band is weaker and bro- 
ader (see Fig. As a consequence, for wave vectors near 
(0, tt) and (tt, 0) this maximum is lost to the foot of the 
spin-polaron peak on crossing the Fermi level. Since the 
spin-polaron band lies somewhat above the Fermi level in 
this region of the Brillouin zone, the situation looks like a 
part of the Fermi surface disappears here and a gap opens 
between the hole energy band and the Fermi level [51l30| . 
In its size (~ 20 meV) and symmetry this gap is similar 
to the pseudogap observed in photoemission spectra. The 
gap disappears at x « 0.12 when the whole Fermi surface 
stems from the crossing of the spin-polaron band with the 
Fermi level. 

The normalized hole density of states, 
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Fig. 5. The normalized hole density of states. 



is shown in Fig. [S] for three hole concentrations from the 
extremely low to the moderate doping case. As seen from 
the figure, a considerable part of the density of states 
is concentrated at a maximum which persists near the 
Fermi level in a wide range of hole concentrations. This 
maximum is produced by the above-mentioned extended 
van Hove singularities in the spin-polaron band and plays 
a great role in the superconductive transition in the t-J 
model |3T| . 



4 The spectrum of spin excitations 

Typical shapes of the spin spectral function B(\&d) are 
shown in Fig. For low and moderate x the spectrum 
consists of a maximum with an extended high-frequency 
tail. Some fine structure can be observed for certain wave 
vectors which is connected with the frequency dependence 
of the polarization operator JBJ. 

The increases of the temperature and the hole con- 
centration act in a similar manner on the spectrum. This 
action is different in the nearest vicinity of (it, tt) and in 
the remainder of the Brillouin zone. As seen from Fig. |6] 
in the latter case with increasing x (and T) the maxi- 
mum in B(\iLu) is shifted to lower frequencies and loses 
its intensity, while in the former case the frequency of the 
maximum, on the contrary, grows. This frequency growth 
is connected with the gap in the spin excitation spectrum 
at (tt, tt). The magnitude of this gap grows with increasing 
x and T. 

As follows from Eq. J5J, the frequencies of spin excita- 
tions satisfy the equation 



ReTT(kcj) - u% = 0. 



(7) 



p(lj) = n- 1 J2M^), 



Their dispersion along the symmetry lines is shown in 
Fig-0for different x and T. In this figure vertical bars de- 
pict the decay widths |Imi7(ko;)|/(2tJk) of the excitations, 
the imaginary parts of the respective poles in Green's func- 
tion. As seen from Fig. Ek, for low x and T the dispersion 
of spin excitations is close to the dispersion of spin waves. 
The main difference is the above-mentioned spin gap at 
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Fig. 6. The spin spectral functions for k = (0, 7r) and (tt, tt) 




Fig. 7. The dispersion of spin excitations. Vertical bars show 
decay widths \lmll (ku))\ / (2^^) . The notations of the symme- 
try points are the same as in Fig. [3] 



(tt, tt). In an infinite crystal the magnitude of this gap is 
directly connected with the spin correlation length £. In- 
deed, using Eq. J5J) and taking into account that the region 
near (tt, tt) gives the main contribution to the summation 
over k, we find for large distances and low temperatures 

/>OC 

(*f*5) = N- 1 £ e* kl / du coth (— ) B(kco) 

oc e^^/lll) 1 ^-^, (8) 

where Q = (tt, tt) and a is the intersite distance (in the 
considered temperature and hole concentration ranges loq 
3> |ReI7(Q, wq)|; therefore in the above equation the gap 
magnitude is approximated by loq). For low x we found 
A 0.2x and consequently 

This relation between the spin correlation length and the 
hole concentration has been experimentally observed in 
La 2 -xSr x Cu0 4 QI|. 

With growing x the spin excitation branch is destroyed 
in some region around the r point - for such momenta 
Eq. (J7J) has no real solution due to a negative value of 
Rc n(kuj). In Fig.0tb.is peculiarity is reflected in the rup- 
ture of the dispersion branch in this region where excita- 
tion frequencies are purely imaginary. Thus, much like the 
Heisenberg model |32U33| . properties of elementary exci- 
tations near (0,0) and (7r,7r) arc different: in the former 
region the excitations are overdamped, while in the latter 
they are gapped. With rise of temperature the branch are 
partly restored around (0, 0). The explanation for such be- 
havior follows from Fig. |H1 For low temperatures and long 
wavelengths Re7T(kw) has a pronounced dip at low fre- 
quencies. This dip is connected with the spin-polaron band 
and is responsible for the lack of real solutions in Eq. Q . 
With increasing T the spin-polaron peaks are smeared, 
the depth of the dip in Fig.|H|becomes smaller and Eq. (JJJ 
has again real solutions. 



5 Magnetic properties 

As already mentioned, in accord with the Mermin- Wagner 
theorem the considered 2D system is in the param- 
agnetic state for T > 0. This result can also be obtained 
using equations of Sec. II. Numerical calculations and the 
analysis of experimental data presented strong evidence 
that the 2D nearest-neighbor s = ^ Heisenberg antiferro- 
magnet has long-range order at T = [SHI ■ Let us consider 
an infinite crystal at T = and find the hole concentration 
x c which destroys this ordering. Notice that the indication 
of the long-range order is a finite value of the condensa- 
tion parameter C which is determined from constraint P| 
and the condition A = 0, 



c =V i 



\Cx\ 1 



' N ^ 



1 - 7k 

l + 7k : 



N 



CO. 



(9) 



8 



A. Sherman, M. Schreiber: Two-dimensional t-J model at moderate doping 



-0.004 



-0.012 



1 I 1 I 1 I 1 I 1 




T"=0.01f 




T 0.0'lf 




r=o.if 




x=0.12 




k=(0,0.1n) 



-u.u io -| 1 1 1 1 1 1 1 1 1 1 

0.0 0.1 0.2 0.3 0.4 0.5 

mlt 

Fig. 8. The real part of the polarization operator Re77(ka;) 
for k = (0, 0.l7r), x — 0.12 and various temperatures. 

As seen from Eqs. (0} and JjJJ, 

C = J^A lim ' 

V a n^oo NVA 

As follows from Eq. (|5J), the sublattice magnetization M = 
(lirri|i|_ ( . 00 |(si • Sq)!) 1 / 2 i s connected with the condensation 
parameter by the relation M = (3C/2) 1 / 2 . To estimate x c 
the dependencies of C\ and a on x have to be taken into 
account. From our calculations for low x we found that 
Id] « 0.2117 - 0.5750a; and a w 1.802 - 8.021x. Substi- 
tuting these values into Eq. © we find that C vanishes at 
x c ~ 0.02. Thus, at T = the infinite 2D crystal is in the 
state with long-range antiferromagnetic order for x < x c 
and in the paramagnetic state for x > x c . 

In the considered finite lattice the correlation length 
is limited by the size of the lattice. The spin correlations 
= (sfso); 1 — (m>n>) calculated from the obtained 
spin spectral function with the use of Eq. JSJl are shown 
in Fig. |21 For large enough x and T the correlations de- 
cay exponentially with distance in the considered finite 
lattice. As mentioned, the method used has no preset 
magnetic ordering. The character of the ordering is de- 
termined in the course of the self-consistent calculations. 
The magnetic susceptibility obtained in these calculations 
is strongly peaked at the antiferromagnetic wave vector 
for low frequencies. No indications of incommensurability, 
which can be related to stripes or other types of phase 
separation, are observed in the susceptibility. Conceivably 
such phase separations are not connected with the strong 
electron correlations described by the t-J model. 

The magnetic susceptibility is connected with the spin 
Green's function J2| by the relation 

X 2 (M = -4 M |D(kw), 

where /is is the Bohr magneton. Experiments on inelas- 
tic neutron scattering give information on the susceptibil- 
ity which can be directly compared with the calculated 
results. Such comparison with the results measured in 
normal-state YZ^CuaOy-j, PI] is carried out in Fig. 1101 




0.00 0.04 0.08 0.12 0.16 

x 




2 4 6 8 10 



m 

Fig. 9. (a) Spin correlations vs. x for T — 0.02t. (b) Spin 
correlations along the diagonal of the crystal [i.e. 1 = (m, m)] 
for x = 0.12. The respective temperatures are indicated near 
the curves. 
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Fig. 10. The imaginary part of the spin susceptibility for k = 
(7r,7r). Curves show calculated results for T = 0.02i ss 116 K, 
x — 0.043 (a) and 0.08 (b). Squares are experimental results 
obtained in normal-state YBa2Cu307_j / at T — 100 K for y — 
0.5 (a) and 0.17 (b) 
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Fig. 11. The local spin susceptibility. The experimental re- 
sults |34| in YBa2Cu30e.5 for uj — 25 meV are shown by filled 
squares. Our calculated data for x — 0.05 and ui — 22.5 meV 
are displayed by open squares. 



YBa2Cu307_j, is a bilayer crystal and the symmetry al- 
lows one to divide the susceptibility into odd and even 
parts. For the antiferromagnetic intrabilayer coupling the 
odd part can be compared with the calculated results. The 
oxygen deficiencies y = 0.5 and 0.17 in the experimen- 
tal data in Fig. ^| correspond to the hole concentrations 
x ~ 0.05 and 0.11, respectively (35]. As seen from Fig. 1101 
the calculated data reproduce correctly the frequency de- 
pendence of the susceptibility, the values of the frequency 
for which Imx(7r, 7r) reaches maximum and their evolution 
with doping. The increase of the frequency of the maxi- 
mum with x reflects the respective growth of the spin gap 
(see Fig. |SJ). In absolute units the calculated maxima of 
Imx(7r, 7r) are 1.5 — 2 times larger than the experimental 
values which is connected with some difference in decay 
widths of spin excitations. 

The temperature variations of the experimental and 
calculated susceptibilities are compared in Fig. ^2 This 
figure demonstrates the imaginary part of the local spin 
susceptibility which is defined as 

Imxi(w) = iV _1 ^lmx(kw). 

k 

As seen from the figure, the calculated temperature vari- 
ation of the susceptibility is also in good agreement with 
experiment. 

The calculated temperature and concentration depen- 
dencies of the uniform static spin susceptibility, 



Xo = x(k^0,o; = 0)=4/i|r- 



are shown in Fig.^J The values lie in the range 2-2.6 eV _1 
which is close to the values 1.9-2.6 eV _1 obtained for 
YBa2Cu307_ y Close values of xo were also obtained 
by the exact diagonalization of small clusters [H]- The 
temperature dependence of xo has a maximum and the 
temperature of the maximum T m grows with decreasing 
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Fig. 12. The uniform static spin susceptibility vs. tempera- 
ture (a) and hole concentration for T — 0.02i (b). 



x. Analogous behavior is observed in cuprates for large 
enough x 37 ,38 ,39 , 40 . In Fig. 112b, T m w 600 K which is 
close to the value observed in L^-zSr^CuO^ for compa- 
rable hole concentrations |37j . As known, in the undoped 
antiferromagnet T m w J 41 j. On the high-temperature 
side Xo(T) tends to the Curie- Weiss behavior 1/T. 

The decrease of \o below T m is sometimes considered 
as the manifestation of the spin gap. In our opinion this 
statement is incorrect. For moderate x and T the long- 
wavelength part of the spin excitation spectrum does not 
feel the gap at (ir,ir). For small but finite values of k 
X(k, 0) oc f™^ du'B(kuj')/uj'. As indicated in Sec. IV, the 
function S(kw') has a maximum which is shifted to lower 
frequencies and loses its intensity with increasing temper- 
ature for such wave vectors. In the above integral the max- 
imum is superimposed with the decreasing function 1/ui' 
which finally leads to the nonmonotonic behavior of \o (T) ■ 

As seen from Fig.ll2h. the two curves for the different x 
are very close in shape and can be superposed by scaling to 
the same values of the maximum xo and T m . Analogous 
scaling was observed in La2_ 2: Sr 2: Cu04 |HZj- As follows 
from the above discussion, the source of this scaling is 
that holes and temperature fluctuations lead in a similar 
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manner to the softening of the maximum in B{ku>') for 
long wavelengths. 

As in experiment 37,38 , the dependence Xo{x) in 
Fig.lT^b has a maximum. However, the calculated value of 
the hole concentration which corresponds to the maximum 
is much smaller than the experimental value x f=a 0.25. A 
possible reason for this discrepancy is the approximation 
used for calculating the hole self-energy 

The spin-lattice relaxation and spin-echo decay rates 
were calculated with the use of the equations 021 



1 



TnI^ ^ (k) — z — ' 



LJ -> 0, 



1 



(63T 2G )2 128 M | N 



[)m 'i^ 63 F e 2 (k)[Re X (kO)] 2 (10) 



where the form factors are 



l^ 63 F e (k)Re X (kO) 



(>:s 



!F|,(k) = (Aj.+4B7 k ) 



(>:i 



! F e (k) - (A,, +4B 7k )' 



(11) 



63 F ± (k) = -[ 63 Fn(k) + 63 F e (k)], 



7 F||(k) = 2C 2 (l + 7k ). 

In the above formulas the hyperfine coupling constants 
B = 3.82 • 10- 7 eV, A ± = 0.845, A\\ = -45, and C = 
0.91B The superscripts a = 63 or 17 indicate that 
the respective quantity belongs to the Cu or O site, re- 
spectively. The subscripts || and _L refer to the direction 
of the applied static magnetic field H with respect to the 
axis c perpendicular to the Cu-0 plane. The form factor 
63 F e is the filter for the Cu spin-echo decay time 63 T2G- 
Due to the different momentum dependencies of the form 
factors ()ll(l measurements of the spin-lattice and spin- 
echo decay rates allow one to extract the information on 
the low-frequency susceptibility in different regions of the 
Brillouin zone. 

Our calculated results are compared with the respec- 
tive experimental data in Fig. 1131 The calculations re- 
produce satisfactorily the main peculiarities of the tem- 
perature dependencies of the spin-lattice relaxation and 
spin-echo decay rates. The growth of ( 63 Ti T) _1 with de- 
creasing x is connected with the increase of the spectral 
intensity of spin excitations near (tt, it) which make the 
main contribution to this rate. For the same hole concen- 
tration ( 63 TiT) _1 is one order of magnitude larger than 
( 17 TiT) _1 . This is a consequence of the fact that Im^ 
is strongly peaked near {%, 7r) and the momentum depen- 
dencies of the form factors l|ll|l . The calculated spin- 
lattice relaxation rates are smaller than the experimental 
values due to the approximation made in the calculation 
of D(kw) which somewhat underestimates Imx at low fre- 
quencies. 




Fig. 13. The temperature dependencies of the spin-lattice 
relaxation and spin-echo decay rates at Cu (a-c) and O (d) 
sites. Open circles with right axes represent experimental re- 
sults, filled circles with left axes are our calculations. (a,c,d) 
Calculations for H|]c and x — 0.12, measurements 139143) in 
YBa2Cu306.63 {x ~ 0.1 |35p. (b) Calculations for nonori- 
ented configuration with x = 0.043, measurements |44) in 
Lai.96Sro.o4Cu04- 



For moderate x with increasing T the low-frequency 
region of Im^(k « Q), Q = (71", 7t) first grows due to 
the temperature broadening of the maximum in its fre- 
quency dependence and then decreases due to the tem- 
perature growth of the spin gap (see Fig. and the re- 
lated discussion) . This nonmonotonic behavior of the sus- 
ceptibility shows up in the spin-lattice relaxation rate at 
Cu in Fig. IT^k. The temperature variations of Im% and 
(^T^T) -1 can be related with the temperature behav- 
ior of the magnetic correlation length. As can be seen in 
Fig. 1141 for moderate x and low T the magnitude of the 
spin gap is determined by the hole concentration and does 
not depend on T. This temperature range corresponds to 
the growth stage in Fig. 113b . As follows from Eq. JSJ), 
the independence of the gap from T means that in this 
temperature range £ does not depend on temperature ei- 
ther which is a distinctive feature of the quantum disor- 
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Fig. 14. The temperature dependency of \f~A oc £ _1 for x = 
0.12. 

dered regime |42II45| . For temperatures above the maxi- 
mum of ( 63 Ti T) -1 we found that 63 T 1 T/ 63 T 2G « const 
(see Figs. 113b and c). As seen in Fig. in this tem- 
perature range £ _1 oc \f~A varies linearly with T . Both 
these facts indicate that in the mentioned temperature 
range the crystal is in the quantum critical z = 1 regime. 
These conclusions are consistent with the phenomenologi- 
cal treatment of experiment in YBa2Cu30j_j, carried out 
in Ref. 05] • The temperature of the maximum in Fig. 113b 
is close to the parameter T* of that work, the tempera- 
ture which separates the quantum disordered and quan- 
tum critical z — 1 regimes. 

For small hole concentrations the temperature range 
in which holes determine the correlation length is ab- 
sent or very small. In these conditions the spin gap grows 
starting from low temperatures and ( 63 TiT) _1 decreases 
monotonously, as shown in Fig. 113b . 

Due to the form factor 17 Fj|(k), Eq. (|ll|l . the momen- 
tum region near (7r,7r) does not contribute to ( 17 Ti T) -1 . 
As indicated in Sec. IV, there is a cardinal difference be- 
tween the behavior of Im \ for k~Q and away from 
(7r, 7r). Due to the spin gap in the former case the frequency 
of the maximum in Imx(w) increases with temperature, 
while in the latter case it decreases. This frequency soft- 
ening leads to the growth of the low- frequency Im \ and 
( 17 TiT) _1 at low T and their saturation for higher tem- 
peratures. Analogous behavior is observed in experiment 
|4*3] . as seen in Fig.fHfli. Thus, the temperature and con- 
centration variations of the spin excitation spectrum in 
the t-J model put forward the simple explanations for the 
behavior of the spin-lattice relaxation and spin-echo decay 
rates observed in cuprates. 



6 Concluding remarks 

In this work we applied Mori's projection operator tech- 
nique and the decoupling of the many-particle Green's 



functions for obtaining the closed set of self-energy equa- 
tions which describes energy and magnetic properties of 
the t-J model of the Cu-0 planes in perovskite high-T c su- 
perconductors. The equations retain the rotation symme- 
try of spin components and zero site magnetization in the 
paramagnetic state which is set for T > or x > x c « 0.02 
in the case T = in an infinite crystal. 

For the parameters of cuprates in the cases of low and 
moderate doping this self-consistently calculated solution 
of the equation is homogeneous. This result indicates that 
phase separations are not related to strong electron cor- 
relations described by the model. 

A number of unusual spectral and magnetic properties 
of cuprate perovskites is satisfactorily reproduced by the 
calculations. Among these properties are the extended van 
Hove singularities around (0, tt) and (tt, 0). Due to strong 
electron correlations these singularities persist near the 
Fermi level in a wide range of hole concentrations which 
has to play an essential role in the superconductivity of 
cuprates. 

The pseudogap in the calculated hole spectrum has the 
same symmetry and is close in magnitude to the pseudo- 
gap observed in photoemission of these crystals. As in ex- 
periment, the calculated pseudogap disappears when the 
hole concentration approaches the optimal doping. The 
existence of two maxima with similar dispersions in the 
photoemission spectrum of lightly doped Ca2Cu02Ci2 is 
related to hole vibronic states in the region of the dis- 
turbed antiferromagnetic order. 

In the considered model the concentration dependence 
of the magnetic correlation length is the same as that ob- 
served in La2- x Sr x Cu04. As indicated, in contrast to this 
crystal, the model does not show low-frequency incom- 
mensurate spin fluctuations for the considered parame- 
ters. 

In other respects the calculated magnetic susceptibility 
is close to that derived from experiments. The results of 
the calculations offer explanations for the observed scaling 
of the static uniform susceptibility and for the changes 
in the spin-lattice relaxation and spin-echo decay rates 
in terms of the temperature and doping variations in the 
spin excitation spectrum. The scaling is related to the fact 
that holes and temperature fluctuations lead in a similar 
manner to the softening of the maxima in the spin spectral 
function for long wavelengths. This softening leads to the 
monotonic increase of the spin-lattice relaxation rate at 
O. 

Contrastingly, due to the spin gap the frequencies of 
the maxima in the spin spectral function for momenta near 
(tt, tt) grow with temperature and hole concentration. This 
region of momenta makes the main contribution to the 
spin-lattice relaxation rate at Cu. For low hole concentra- 
tions the growth of the frequencies of the maxima leads to 
the monotonic decrease of the rate. For moderate concen- 
trations two temperature regions may be distinguished. 
In the low-temperature region the gap magnitude and the 
correlation length depend only weakly on temperature. 
In this region the spin-lattice relaxation rate at Cu grows. 
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The behavior of these parameters in the high-temperature 
region is the same as for low hole concentrations. 

This work was partially supported by the ESF grant 
No. 4022 and by DFG. 
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